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Abstract 

We develop a mathematical theory to represent dislocations and discli- 
nations in single crystals at the mesoscopic scale by considering con- 
centrated efTects, governed by the distribution theory, combined with 
multiple-valued kinematic fields. Our approach provides a new under- 
standing of the continuum theory of defects as developed by Kroner (1980) 
and other authors. The fundamental identity relating the incompatibility 
tensor to the Frank and Burgers vectors is proved in the 2D case under 
appropriate assumptions on the strain curl behaviour in the vicinity of the 
assumed isolated defect lines. In general our theory provides a rigorous 
framework for the treatment of crystal line defects at mesoscopic scale. 
Eventually this work will represent a basis to strengthen the mathematical 
theory of homogenization from mesoscopic to macroscopic scale. 

1 Introduction 

Dislocations can be considered as the most complex class of defects for several kinds 
of single crystals (Volkl & Miiller 1994; Jordan et al. 2000) and the development of a 
relevant and accurate physical model represents a key issue with a view to reducing the 
dislocation density in the crystal by acting on the temperature field and the solid-liquid 
interface shape during the growth process (Dupret and Van den Bogaert 1994). 

However the dislocation models available in the literature, such as the model of 
Alexander and Haasen (1986), are often based on a rather crude extension of models 
initially developed for polycrystals (as usual metals and ceramics are). In this case. 
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some particular features of single crystals, such as material anisotropy or the existence 
of preferential glide planes, can be taken into account up to some extent, but the 
fundamental physics of dislocations in single crystals cannot be captured. In fact, 
dislocations are lines that either form loops, or end at the single crystal boundary, or 
join together at some locations, while each dislocation segment has a constant Burgers 
vector which exhibits additive properties at dislocation junctions. These properties 
play a fundamental role in the modelling of line defects in single crystals and induce key 
conservation laws at the macro-scale (typically defined by the crystal diameter). On 
the contrary, no dislocation conservation law exists at the macro-scale for polycrystals 
since dislocations can abruptly end at grain boundaries inside the medium without 
any conservation law holding across these interfaces. 

Aware of these principles and of the pioneer works of Volterra (1907) and Cosserat 
(1909), Burgers (1939), Eshelby (1956, 1966), Eshelby, Frank & Nabarro (1951), Kondo 
(1952), Nye (1953), and Kroner (1980) among other authors (Bilby 1960; Nabarro 
1967; Mura 1987) consider a tensorial density to model dislocations in single crystals at 
the macro-scale, in order to take into account both the dislocation orientation and the 
associated Burgers vector (cf. the survey contributions of Kroner 1980, 1990, Kleinert 
1989 and Maugin 2003). However, in these works, the relationship between the macro- 
scale crystal properties and the basic physics governing the nano-scale (defined by 
the inter-atomic distance) is not completely justified from a mathematical viewpoint. 
Therefore, to well define the concept of tensorial dislocation density, we here introduce 
the meso-scale as defined by some average distance between the dislocations. The 
laws governing the dislocation behaviour are modelled at the nano-scale, while the 
meso-scale (defined from the nano-scale by ensemble averaging or by averaging over 
a representative volume (Kroner 2001)) defines the "dislocated continuous medium", 
where each dislocation is viewed as a line and the interactions between dislocations 
can be modelled while the laws of linear elasticity govern the adjacent medium. 

The present paper focuses on meso-scale modelling with a further view to clarifying 
the homogenization process from meso- to macro-scale. This latter issue is addressed 
in the companion work of Van Goethem & Dupret (2009a). Since dislocations are 
lines at the meso-scale, concentrated effects must be introduced in the mesoscopic 
model as governed by the distribution theory (Schwartz 1957). In addition, since 
integration around the dislocations generates a multiple- valued displacement field with 
the dislocations as branching lines, multivalued functions must be considered (cf, e.g., 
Almgren 1986). This combination of distributional effects and multivaluedness is a key 
feature of the dislocation theory at the meso-scale but unfortunately the difficulties 
resulting from this mathematical association have not well been addressed so far in 
the literature (see also Thom 1980). As an example, non-commuting differentiation 
operators are freely introduced without any justification by Kleinert (1989). Therefore, 
the principal objective of this paper is to provide a strong mathematical foundation 
to the meso-scale theory of dislocations, showing how the distribution and geometric 
measure theories can be correctly used with multiple- valued fields. 

In fact, a key modelling issue arises from the fact that homogenization from meso- 
to macro-scale has no meaning for multiple-valued fields such as displacement and 
rotation, since this operation is exclusively allowed for additive (or extensive) fields 
such as stress, energy density or heat flux. This observation becomes obvious when 
homogenization is defined by an ensemble averaging procedure, since multiple- valued 
fields are mathematically defined as extended functions which cannot be added since 
their "domains" depend on the defect line locations. This consideration justifies the 
present analysis. For the sake of generality, disclinations, which represent a second 
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but rarer kind of line defect, with in addition a multiple- valued rotation field, are here 
considered together with dislocations. 

In the literature the macroscopic dislocation density is classically defined as the 
curl of the plastic distortion (Ifead et al. 1993; Cermelli & Gurtin 2001; Koslowski 
et al. 2002), following a postulated distortion decomposition into elastic and plastic 
parts. However, this decomposition cannot be rigorously justified (contrarily to the 
strain decomposition) since elastic and plastic rotations cannot be set apart without 
some hidden arbitrariness. In contrast, the present paper paves the way for a rigourous 
definition and treatment of the macroscopic dislocation density, as obtained from well- 
defined mesoscopic fields under precise geometric-measure model assumptions, and 
from which the distortion decomposition can be obtained together with its relationship 
with the dislocation density (Van Goethem & Dupret 2009a). 

The present paper is restricted to the mesoscopic 2D theory for a set of assumed 
isolated dislocations and/or disclinations. This theory is extended to the case of count- 
ably many dislocations in Van Goethem & Dupret (20096) where the appropriate 
mathematical objects and functional spaces are ultimately defined for homogenization 
to the macro-scale. This latter paper will be referred to as Part B in the sequel. Ex- 
tension to the dynamic 3-D case is under investigation. Eventually, the complete link 
between the mesoscopic and macroscopic behaviours of single crystals with line defects 
should be derived from these developments. In tlie scaling analysis summarized in 
this introduction is detailed and the basic concepts used to represent the dislocated 
continuous medium are introduced. The general mathematical theory is developped 
in 331 while in the 2D distributional theory of the dislocated continuous medium is 
established in the case of isolated parallel dislocations/disclinations. Conclusions are 
drawn in ^ 

2 Multiscale analysis of dislocations 

2.1 Nano-scale analysis: crystalline lattice 

At the nano-scale the characteristic length is the interatomic distance and the reference 
body is a perfect lattice. Given a dislocation in the general sense (dislocation and/or 
disclination) , the atomic arrangement at time t generaly differs from the reference 
arrangement, but however the atom displacements are not uniquely defined (Kleinert 
1989). Indeed any atom of the reference configuration can in principle be selected to 
define the displacement of a given atom of the actual configuration which therefore is 
a multivalued discrete mapping. Moreover, in general, the dislocation position cannot 
be determined precisely at the atomic level since several dislocation locations in the 
actual crystal can be associated with the same picture of the atom positions. In fact 
the defect should be understood as located inside a nanoscopic lattice region. 

2.2 Meso-scale analysis: dislocated continuous medium an 
associated reference configurations 

At the meso-scale the characteristic length is some average distance between two neigh- 
bour dislocation lines. This scale is the one on which this paper focuses, in the frame- 
work of 2D linear elasticity. At time t, the body is referred to as TZ*{t) as corresponding 
a random sample corresponding to a given growth experiment. 
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A reference configuration TZq with respect to the actual configuration 7l*{t) is any 
selected one-to-one transformation of TZ*{t). TZg may be chosen as being the body 
at any given (past or future) time to or in contrast be a fictitious transformation of 
7l*{t), and the displacement and rotation fields on TZ*{t) {u* and ul) are then defined 
with respect to the chosen TZq. In the present TZq will always be defined as stress-free 
and without dislocations, but its selection will remain arbitrary up to this restriction 
and hence (and this is a keypoint) the defect governing laws must be invariant with 
respect to the choice of TZq . 

It will be precised later that displacement and rotation are multivalued fields at 
the mesoscale, and hence are defined on a set called a Riemann foliation F (and not 
of T^-o). The set F can be univoquely associated to TZ*{t) if a cut is introduced in the 
foliation in order to select one particular branch of the displacement and rotation. 

In view of multivaluedness and the existence of a family of acceptable reference 
configurations, the main field of our analysis is the assumed linear elastic strain which 
is clearly single valued and independent of the choice of TZq. The Burgers vector B* 
and Frank vector fl* are key invariant quantities related to the jump of the multival- 
ued displacement and rotation fields, as directly derived from the linear strain. Their 
precise definition will be given in i|3] 

At this stage, some definitions and assumptions have to be introduced. 

Notations 2.1 In the following sections, the assumed open domain is denoted by Q, 
(in practice but not necessarily Q is bounded), the defect line(s) are indicated by C dO,, 
and flc is the chosen symbol for which is also assumed to be open. 

Definition 2.1 (3-D mesoscopic defect lines) At the meso-scale, a 3D set L of 
dislocations and/or disclinations is defined as a set of isolated rectifiable arcs L^''\ k £ 
X, without multiple points except possibly their extremities and on which the linear 
elastic strain is singular. Here a set of isolated arcs means a set of arcs: (i) whose 
extremeties form a set of isolated points of fl in the classical sense and (ii) such that 
each point x of these arcs except their extremities can be located in a smooth surface 
S{x) bounded by a loop C{x) and such that S{x) \ x € Qc ■ 

Assumption 2.1 (3D mesoscopic elastic strain) Henceforth we will assume that 
the linear strain is a given symmetric L^{Q) tensor ij prolonged by on the 

dislocation set L and compatible on Q.c ■ In other words, the incompatibility tensor, as 
defined by 



where derivation is intended in the distribution sense, is assumed to vanish everywhere 
on Qc- 

Let us now introduce the dislocation and disclination density tensors (A*j and 0*^) 
which are the basic physical tools that will be used to model defect density at the 
meso-scale. 

Definition 2.2 (defect densities) 




(2.1) 



DISCLINATION DENSITY: 




(2.2) 



feeicN 



DISLOCATION DENSITY: 




(2.3) 



fceicN 



'^It should be noted that L'J(Q) with any I < q < 2 would hold as well. 
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where symbol Sj^{h) is used to represent the one- dimensional Hausdorff measure^ den- 
sity concentrated on the rectifiable arc L'*'-' with the tangent vector r^''^ defined almost 
everywhere on L^^\ while n**''^ and B*^''^ denote the Frank and Burgers vectors of 
L^^\ respectively. 

The present paper and Part B address the 2D problem only. Generalization of our 
theory to the 3-D case will be considered in further publications. 

Definition 2.3 {2D mesoscopic defect lines) At the meso-scale, a 2D set C of 
dislocations and/or disclinations is defined as a set of isolated parallel lines L^^\ i £l, 
on which the linear elastic strain is singular. In the sequel, these lines will be assumed 
as parallel to the z-axis. 

More complex sets of 2D defect lines are considered in Part B. 

Remark 2.1 The term 2D here refers to the structure of the countable union of 
points, denoted by lo, located at the intersection between £ and the z = zg-plane. 
In this context, the strain is said 2D if it solely depends upon the coordinates Xa £ 
(a = 1,2). In that case, the displacement and rotation fields will generaly depend on 
the three space variables. 

Notations 2.2 (2D defect densities and incompatibility) In2D, the vectors r]l,Q*^. 
and K% will denote the tensor components rj*^, QJ^. and K*f.. Greek indices will be used 
to denote the values 1,2 (instead of the Latin indices used in 3D to denote the values 
1,2 or 3). Moreover, eajs will denote the permutation symbol e^ap. 

The disclination and dislocation density tensors Q^. and will be shown in this paper 
to be related by a fundamental distributional relation to the strain incompatibility 77^. 
In fact, under suitable assumptions on the strain curl (the so-called Frank tensor), the 
following theorem will be proved in the 2D linear elastic case. 

Main theorem: 

incompatibility decomposition for 2D isolated defect lines. 

The mesoscopic strain incompatibility for a set of isolated parallel rectilinear dislo- 
cations L writes as 

nl = 'dl + eo,i3danlp, (2.4) 
where denotes the contortion tensor, 

i^li3 = Skzap -^OL*,5ki3, (2.5) 
with standing for an auxiliary defect density vector, 

al := Al - Skaea/jQlixii - xois), (2.6) 
and where xo is a selected reference point in Q,. 



■^The reader is referred to Mattila (1995) for details on HausdorfT measures. 
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The latter result appears in Kroner's work (1981) under assumptions which are 
not compatible with our approach. In fact, in his work this result follows in a straight- 
forward manner from an "elastic-plastic" displacement gradient (or distortion) decom- 
position postulate, which itself requires the selection of a particular reference config- 
uration and does not properly handle the intrinsic multivaluedness of the mesoscopic 
problem. Moreover, in our result the link between the defect densities and the Frank 
and Burgers vectors is clearly made, and precise assumptions on the strain field and 
the admissible defect structures are provided in order to validate the result. 

The above theorem will be generalized in Part B to the case of a countable union of 
parallel rectilinear dislocations. Eventually, the required "single- valued" distributional 
fields will be defined in the appropriate functional spaces for their homogenization to 
the macro-scale. 

3 Multiple- valued fields and line invariants in 
the 3D case; distributions as a modelling tool 
at the meso-scale 

Notations 3.1 In the following sections, x or Xi will denote a generic point of the 
defect line(s), x or Xi a generic point ofQ,c, and xq or xoi a given fixed reference point 
of Q,c ■ When x and x are used together, x denotes the projection of x onto a given 
defect line in an appropriate sense and Vi :— Ui{x,x) is the unit vector joining x to 
x. The symbol 0e is intended for a set of diameter 2t enclosing the region C. More 
precisely, 0e is defined as the intersection with Q. of the union of all closed spheres of 
radius e centred on C: 

0£ := n n [J B[x,e]. 

If C consists of an single line L, 0^ is a tube of radius e enclosing L. 

Notations 3.2 In the sequel, considering a surface S of Q crossed by a dislocation 
L at X and bounded by the curve C , symbols dC , dL, and dS will denote the ID 
Hausdorff measures on C and L, and the 2D Hausdorff measure on S, respectively, 
with Tj standing for the unit tangent vector to L at x (when it exists). In some cases 
(having fractal curves in mind) the symbols dxk and dSi := eijkdx'^p dx^^^ will stand for 
infinitesimal vectors oriented along C and normal to S, respectively, with m addition 
dCi{x) := (.imndXmTn denoting an infinitesimal vector normal to C when t„ — t„(x) 
exists. 

In the present section, the strain is assumed to satisfy assumption 12.11 and to be 
smooth away from jC. 

3.1 Distributional analysis of 3D multiple-valued fields 

In general, a multivalued function from Q,c to is defined as consisting of a pair of 
single-valued mappings with appropriate properties: 

F^flc and F ^ R^, 
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where F is the associated Riemann foUation (Almgren 1986). In the present case 
of meso-scale elasticity, we will limit ourselves to multivalued functions obtained by 
recursive line integration of single- valued mappings defined on Q.c- Reducing these 
multiple line integrals to simple line integrals, the Riemann foliation shows to be the 
set of equivalence path classes inside Q,c from a given xo £ Q,c with homotopy as 
equivalence relationship. Accordingly, a multivalued function will be called of index n 
on Q,c if its n-th differential is single- valued on Q.c- No other kinds of multifunctions 
are considered in this work, whether £ is a single line L or a more complex set of 
defect lines (with possible branchings, etc.). 

Notations 3.3 The symbol is used for partial derivation of a single- or multiple- 
valued function whose domain is restricted to Qc- Locally around x £ Q,c, for smooth 
functions, the meanings of 9^"' and the classical dj are the same, whereas on the 
entire SI the partial derivation operator dj only applies to single-valued fields and must 
be understood in the distributive sense. A defect-free subset U of is an open set such 
that (7 n £ = 0, m such a way that and dj coincide on U for every single- or 
multiple-valued index-1 function. 

In the following essential definition generalizing the concept of rotation gradient 
to dislocated media, the strain is considered as a distribution on SI. 

Definition 3.1 (Frank tensor) The Frank tensor dmO^t defined as the following 
distribution on S7; 

dmUJif. := ekpqdpSq^, (3.1) 

in such a way that 



< dmOjl,ip >:= ~ I ekpqEqmdptpdV, (3.2) 

J (I 

with a smooth test-function with compact support in Q,. 

In fact, in the vicinity of a defect line the tensorial distribution dmUi'^. is the finite 
part of an integral when acting against test-functions. Indeed, since dp£q^ might 
be non L'^ (r2)-integrable in view of its possibly too strong singularity near the defect 
lines, instead of being directly calculated as an integral, < eupqdpEqm,^ > must be 
calculated on Q, as the limit 

lim / etpqdpSq^ipdV -f / eupqSlm'^dSp , (3.3) 

where the second term inside the parenthesis is precisely added in order to achieve 
convergence. One readily sees after integration by parts that (|3.3p is equal to (|3.2[) 
provided lim f2 \ 0e = Q,c (which is a general hypothesis limiting the acceptable defect 

e— ^0 

lines and certainly holds true for the lines satisfying definition [2TT]) . 

Considering the possibly index-1 multivalued rotation vector ujI, it should be ob- 
served from definition 13.11 that dmij^t = dm^ujk on flc as a consequence of the classi- 
cal relationship between infinitesimal rotation and deformation derivatives. However, 
dmi^k is defined by (|3.1[) as a distribution and therefore concentrated effects on £, and 
its infinitesimal vicinity have to be added to dm^ui^, justifying the use of the symbol 
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dm^k instead of dmUJ^ without giving to dm the meaning of an exact derivation op- 
erator. In particular, it may be observed that the identical vanishing of din'ujl on Qc 
does not necessarily imply that the distribution dm^t vanishes as well. In fact from 
1)3. 3[) . it can be shown in that case that 

< dmiol,ip >^ lim ekpqSqmfdSp = - ekpg£*,ndpipdV, (3.4) 

which is generally non-vanishing. Finally, as soon as the definition of the tensor 
distribution dm^t is given, so are the distributional derivatives of dmO^l'- 

< dldmi^l,(p >= - < d,niul,dlip >^ / ekpnSmndpdlifdV. (3.5) 

Jn 

3.2 3D rotation and displacement vectors 

The rotation vector is defined from the linear strain together with the rotation at a 
given point xq. From this construction follows an invariance property of as a mul- 
tifunction (recalling that multivaluedness takes its origin from the existence of defect 
lines which render the strain incompatible on the entire fi). 



Starting from the distributive definition 13.11 of dmi^l^ the form dm^^k d^m is in- 
tegrated along a regular parametric curve F C Qc with endpoints Xq,x G Qc- For 
selected xo and cjofei the multivalued rotation vector is defined as0 



d^m, 

where #F is the equivalence class of all regular curves homotopic to F in Q,c- Indeed, 
from strain compatibility in Q.c, i.e. from relation (|2.ip . it is clear that oj^ is a function 
of #F only. Consider now a regular parametric loop C (in case C is a planar loop, it 
is a Jordan curve) and the equivalence class #C of all regular loops homotopic to C 
in Q.C- Here, the extremity points play no role anymore and two loops are equivalent 
if and only if they can be continuously transformed into each other in Q.c- The jump 
of the rotation vector lj^ along #C depends on #C only and is calculated as 

K] = [^1]{#C) = I dmi^ld^m- (3.6) 

Jc 

The following developments address the displacement field multivaluedness as a 
mere consequence of strain incompatibility. The procedure defining the displacement 
vector from the rotation vector by means of line integrals is classical in linear elasticity. 
The following tensor plays in the construction of the displacement field a role analogous 
to dm^k in the construction of the rotation field. 

Definition 3.2 (Burgers tensor) For a selected reference point xq G ilc, the Burg- 
ers tensor is defined on the entire domain Q as the distribution 

dibl{x-xo) := £kii^) + <^kpq{xp - xop)diLu*{x). (3.7) 



•^For a non-smooth strain, integration is to be understood in the distribution sense. 

*We note that C could be non rectifiable, i.e. of infinite length. Integrals on fractal curves 
and the related Stokes' and Gauss-Green's theorems are analysed by Harrison & Norton 
(1992), where it is shown by the C°° smoothness of the difi'erential form dm'^k'^Xm on f2£ 
that Il3.6| l still holds even when the Hausdorff dimension of C is higher than 1. 



8 



The Burgers tensor can be integrated in the same way as the Frank tensor along 
any parametric curve F, providing for selected Uq^. and Uqj, the index 2-muItivalued 
displacement vector u^: 

Ul = ul[x,ij^T-,l^l,ul) = ulk+ eklrni^*l{ij^V-,U}l)[Xra - XOra) + j dlbl{(.)d£,l , 

which is a function of x and #r only (this following from (|2.1|) and (|3.7[) V It may be 
observed that dibl and the vector 

K = 6fc(#r;tlo) =ul~ €klmi^t{^rn - Xq^) (3.8) 

are related in the same way as dmi^l ^-nd ujI, including the fact that dibl = dl'^^bX on 
Qc- The jumps of bl along #C and of ttj at x along #C (which depend on #C only) 
are calculated as 

[&fc](#C; xo) = [ul]{x- i^C; xq) - efcimM*](#C)(2;m - xom) = / dibld^i. (3.9) 

Jc 

Let us now focus on the case of a given isolated defect line I/''-*, i £ I. The jump 
[uil] of the rotation vector ujI around L''-* is defined as the jump of ujI along ^C{x), 
with X a point of i'*^ and C{x) a loop enclosing once the defect line L^*^ and no other 
defect line as specified in definition [2TT] It turns out that this jump is the same for any 
X and suitable C{x). Similarly, the jump [6^] of the vector 6J around I/*-'' is defined as 
the jump of bX along #C{x) and is also the same for any x and suitable C{x), given 
Xo- In fact, the following result is well-known (Kleinert 1989). 

Theorem 3.1 (Weingarten's theorem) The rotation vector uj^. is an mdex-1 mul- 
tifunction on 0,c whose jump OjJ :— [ojI] around the isolated defect line i'^', i £ 2, is 
an invariant of this line. Moreover, for a given xo, the vector b^ is a multifunction of 
index 1 on whose jump Bl :— [b^] around L''-* is an invariant of this line. 

From this result, the Frank and Burgers vectors are defined as invariants of L^^K 

Definition 3.3 (Frank and Burgers vectors) The Frank vector of an isolated de- 
fect line L^'' , i £ I, is the invariant 

Ql ~ H], (3.10) 
while for a given reference point Xo its Burgers vector is the invariant 

Bl ■■= [bl] = [ul](x) - eklm^liXm - XOm). (3.11) 

A defect line with non-vanishing Frank vector is called a disclination while a de- 
fect line with non-vanishing Burgers vector is called a dislocation. Clearly a disclina- 
tion should always be considered as a dislocation by appropriate choice of xq while 
the reverse statement is false since 0,% might vanish. In fact, two distinct reference 
points Xo and x'o define two distinct Burgers vectors, obeying the relation B^ — B'^ = 
ekim{xom — x'om,)^* (notiug that B^ill is an invariant independent of the arbitrary 
choice of 2:0). Therefore, for a non-zero Frank vector, the vanishing of the Burgers 
vector depends on the arbitrary choice of xo. This is why in the present paper, the 
word "dislocation" means in the general sense a dislocation and/or a disclination. A 
pure dislocation is a dislocation with vanishing Frank vector. 
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Remark 3.1 It should be emphasized that the assumption of isolated defect lines is 
required to construct appropriate enclosing loops in order to define their Frank and 
Burgers vectors. In Part B this assumption will be removed for a countable set of 
parallel defect lines under appropriate assumptions on the Frank tensor. 

In general, every defect line will contribute to the rotation and displacement mul- 
tivaluedness, and hence these latter fields are defined over ^Ic and do not share the 
structure of a vector space. In other words, the displacement and rotation fields cannot 
be added since their domains depend on the defect line locations. 

Therefore, besides the strain field which is the seminal ingredient of the present 
theory, the Burgers and Frank tensors appear as fundamental quantities able to char- 
acterize the amount of defects on each single line or in the whole dislocated crystal. 
Together with the geometry of the defect set, these vectors provide the key defect 
measures called the dislocation and disclination density tensors (which now belong 
to a vector space). Accordingly, the following well-known result can be readily shown 
and is fundamental in the framework of our investigations since it implies conservation 
laws at the meso- and macro-scales. 

Theorem 3.2 (conservation laws) Single disclination and dislocation lines are al- 
ways closed or end at the boundary of Q. Moreover, in all cases, 

diO*j = (9iA* = 0. 



4 Distributional analysis of incompatibility for 
a single rectilinear dislocation 

4.1 The 2D model for rectilinear dislocations 

2D elasticity means that the strain £*j is independent of the "vertical" coordinate z. 
However this assumption introduces no restriction on the dependence of the multiple- 
valued displacement and rotation fields upon z. 

Notations 4.1 In this 1\ the single defect line L is assumed to be located along 
the z-axis. The two planar coordinates will be denoted by {x,y) or Xa. The projection 
of X = {xa,z) on L is X — (0,0,^) and (uajO) stands for the unit vector from x to 
X. Symbols (6^,6^,62) or {ea,ez) denote the Cartesian base vectors, while (6^,66,6^) 
denote the local cylindrical base vectors. For a planar curve C , the notation dCa{x) = 
tafidxji is used for an infinitesimal vector parallel to the curve normal. 

Let us observe that many fields are singular at the origin and that Ql is in fact 
the domain where the laws of linear elasticity apply. Moreover, the strain can be 
decomposed into three tensors: 




planar strain 3D shear pure vertical compression / dilation 

Lemma 4.1 (2D compatibility) InQ.L, from 2D strain compatibility, there are real 
numbers K, aa and b such that 

ea-/ei3sdadis£*g — 0, 

eafjdaSp,, = K, (4.1) 
£zz = a^Xa + b. 
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This lemma is easily proved from assumption 12.1 



□ 



Remark 4.1 The present theory does not make use of the linear elasticity constitutive 
laws and the momentum and energy conservation laws, since in the framework of 
Continuum Mechanics arbitrary body forces and heat supply can always be applied to 
the medium. Moreover, the sum of these two body contributions and the unsteady terms 
governing the medium dynamics can generally be nonsmooth, and hence the stress and 
heat flux derivatives have to be treated as mathematical distributions thereby providing 
a physical justification to our approach. 

The remaining of this section will be devoted to present the three classical examples 
of 2D line-defects for which the medium is assumed to be steady, body force free and 
isothermal (detail is given in Van Goethem 2007). 

• Pure screw dislocation. The displacement and rotation vectors write as 



Ui Ci = 



2^ 



■Cz and bJiCi 



-V X uld 



47rr 



(4.2) 



in such a way that the jump [u*] vanishes identically, while the Cartesian strain 
is divergence-free on Q. and writes as 



-Bl 
47rr^ 



y 



(4.3) 



Moreover, inside the Frank tensor writes as 



-Bt 



cos 261 sin 261 
sin 29 - cos 29 




Pure 



47rr^ 

dislocation. The displacement vector is 
u*ei 



(4.4) 



g.^gogj + l) ^B19 
2-K " + 27r ^' 



while the rotation uj* vanishes together with its jump. The Cartesian strain 
(which requires additional regular terms to correspond to balanced stresses) 
writes as 



"By 

27rr2 



X 


y 





y 


—X 















(4.5) 



noting that the tensor d-mOJ% vanishes identically inside 
Wedge disclination. The rotation vector is 

ni9 

with the multiple-valued planar displacement field given by 



[$(l-^*)yln(^)-§(l + ^')x- + §x^ 



'In 



ye 



(4.6) 
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(where v* stands for the 2D Poisson coefficient, to be considered as an arbitrary 
constant togetlier with R) and with a vanishing Burgers vector: 



— iB,, 



- i[u*y] + + ix) = 0. 



The Cartesian strain writes as 



and hence 



47r 



Stt 



(log ^ + 1) 

(log^ + 1) 



" cos 26) sin 26* ' 

sin 29 - cos 261 





2nr 









sin 9 
— cos 6 




(4.7) 



(4.8) 



Remark 4.2 The energy density £* — ^a*j£*j is not -integrable for both kinds of 
dislocations, while it is finite for the wedge disclination. Therefore, a Hadamard finite 
part (Schwartz 1957; Estrada & Kanwal 1989) is needed to represent the compliance at 
the meso-scale (another approach makes use of strain mollification by a so-called core 
tensor (Koslowski et al. 2002)). This issue, whose solution requires to develop matched 
asymptotic expansions around the singular line in accordance with the infinitesimal 
displacement hypothesis, will not be addressed further in the present paper which only 
focuses on the geometry of dislocations. 

Remark 4.3 The above expressions of dislocations and disclinations do not necessar- 
ily provide balanced stresses. The present theory is fully independent of any dynamical 
assumption and only focuses on the geometrical concentrated properties of the defect 
lines. 



4.2 Mesoscopic incompatibility for a single defect line 

For 2D problems the incompatibility vector contains all the information provided 
by the general incompatibility tensor. The latter expresses on the one hand the non- 
commutative action of the defect line over the second derivatives of the rotation vector 
and on the other hand is related to concentrated effects of the Frank and Burgers 
vectors along the defect line. 

Definition 4.1 (2D incompatibility tensor) In the 2D case, the mesoscopic in- 
compatibility vector is defined by 

Vk ■= eapdadpul. (4.9) 

A strain field is compatible if the associated incompatibility vector vanishes. 

As shown in the following sections, concentration effects will be represented by 
means of first- and second-order distributions. 
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Notations 4.2 Let, with use of notation \4- 1\ fiz ~ {x £ Q with a prescribed z} and 
$7? :— i^z \L while the radius r = \\x — x\\ is the distance from a point x inside fl to 
L. Then, the ID Hausdorff measure density concentrated on L will be denoted by 5l- 

In what follows the hypothesis consists in assuming that the strain radial depen- 
dence in the vicinity of L is less singular than a critical threshold. This is verified for 
instance by the wedge disclination whose strain radial behaviour is O(lnr) and by 
the screw and edge dislocations whose strains are 0{r^^). 

For a straight defect line L, according to these examples, the hypotheses on the 
strain and Frank tensors read as follows. 

Assumption 4.1 {2D strain for line defects) The strain tensor £*j is indepen- 
dent of the coordinate z, compatible on Q,l = Q\L in the sense that conditions 114 -1^ 
hold, smooth on f2_L, and -integrable on SI. 

Assumption 4.2 (local behaviour) The strain tensor £*j is assumed to be o{r~^) 
(e — !■ 0"'") while the Frank tensor is assumed to be o{r~^){e — >■ 0^). 

The two following lemmas are needed for the proof of our main result for a single 
isolated defect line. 

Lemma 4.2 Let Ce(£), e > 0, denote a family of 2D closed rectifiable curves. Then, 
in 2D elasticity, the Frank tensor and the strain verify the relation 

hm / (xadi3Uj*dxfi + eK/jfL) dx^ = 0, 

provided the length ofCe is uniformly bounded and as long as the convergence Ci{x) — >■ 
X is understood in the Hausdorff sense, i. e. in such a way that 

max{||a:: — £||, s £ Ce(fi)} — > 0. 
Proof. The second compatibility condition of (|4.1|) is equivalent to 

d-/£0z — 9/3£^7z = Ke^/i, 

from which, in 2D elasticity: 
and 

{XadpLol, + SapeK'yE*^) = 9,8 (^Xae^.-^E*^ - XaKS^/j. 

Since, under the assumptions of this lemma, 

lim / XadxK. = 0, 

while the strain is a single-valued tensor, the proof is achieved. □ 

function /(e) is said to be O (e — > 0+) if there exists K, to > s.t. < e < eo ^ 

< K\g{e)\. A function /(e) is said to be o{g(e)) (e 0+) if lim = 0. 
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Lemma 4.3 In 2D elasticity the planar Frank vector vanishes. 
Proof. Since 

9^6* = Epr + er^{x^ — XQ^)5fiU}l — erj{z — Zo)S0iU*, 

the planar Burgers vector simply writes as 

B* = [£0^ + eT^{x^ — xo^)Si3Uj*) dxp — eT^{z — zo)0.*, 
Jc 

where C is any planar loop. By Weingarten's theorems the Burgers vector is a constant 
while the integrand is independent of 2, from which the result obviously follows. □ 



Theorem 4.1 (main result for a single defect line) Under assumptions \4- 1\ and 

for a dislocation located along the z-axis, incompatibility as defined by equation 
\4-!^ is the vectorial first order distribution 



Vk = SkzV*z + Sk^ijl, (4.10) 

with 

ril = VI^Sl + ea-i {B* - e^^xo/ifll) OcSl, (4-11) 

= ie«aB:9«<5z,. (4.12) 



Proof. For some small enough e > and using notations 13.11 a tube 0e can be 
constructed around L and inside SI. Assuming that the smooth 3D test-function (p 
has its compact support containing a part of L, Q^^z denotes the slice of the open 
n \ 0E obtained for a given x £ L, i.e. 

f^E,z {x G ^Iz such that Ijx-dl > e}, 

while the boundary circle of Qi,z is designated by Ce,z- 

▲ Let us firstly treat the left-hand side of 1)4. 10|) . From definitions 14.11 and 13.11 and 
equations (|3.H) . (|3.2|) and p.3p . it follows that 



where 



<r]l,ifi>= / dz lim Tlkiz,if,€), (4.13) 

E^0 + 



Ylk[z,ip,e) ■— — I ea/idpujldafdS - / eapek^n^pnda'fidC^. (4.14) 

From notation 14.11 the boundedness of \dTds'fi\ on Ql provides the following Taylor 
expansions of if and daif around x: 

ifi{x) = (fi{x) + ruadaip(x) + Y'^Ti^sdrdsifix + '-fi{x - x)) , (4.15) 

daifiix) = daip{x) + ruTdrdaifi {x + ■y2{x — x)) , (4.16) 
with < 71 (x — i), 72(3; — i) < 1. 

▲ Consider the first term of the right-hand side of (|4.14|l . noted life. By virtue of the 
strain compatibility on f2_L and Gauss-Green's theorem, this term writes as 

flk{z,ip,e) ■— - [e-fftdpulip) dS = / e^pd jjulipdC^. 
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Since by notations 14.11 and 14.21 rva ~ Xa — Xa = Xa, then equation (|4.15p and 
assumption 14.21 show that, for e — i- 0"*", 



life = y e^pdi^ul (^ip{x) + Xada^pix)^ 



dC-y + o{l). 



▲ Consider the second term of the right-hand side of (|4.14p . noted life. On account of 
assumption 14.21 and from expansion (|4.16[l . this term may be rewritten as 

= -daip{x) / eall€k-,nSp„dC-, + o{l). 

▲ From Weingarten's theorem and recalling that dC-y — e-yrdxr, the expression life — 
life + life then writes as 

life = daifi{x) / [xcOtUjI - eapek-fne^T-San) dxT 

+Ql^{x) + o{l). (4.17) 

▲ Consider the first term of the right-hand side of (|4.17p . noted Hj., and take ^ = 7 
in the identity 

^k^n^lT = Skz (Sj^Snr - S^jS^^) — 5nz {SjfSkr " SkjS^f) (4.18) 

in such a way that 

ll'k =da'p{x)j [Xadr^l - 5kz^al3£lT + 5kT^afi£*fiz) dx^. (4.19) 

▲ The cases k — z and k — n are now treated separately. 

• When k — z, definition 13.21 shows that 



:— Epr + (-T',{x~, — xot)di3Lol — eT^{z — zo)d0U}* 

which, after multiplication by era and using (|4.18|) with r, a and z substituted for k, ^ 
and n, is inserted into (|4.19|l . thence yielding: 

n!, = datp{x) / (eradpb* + xoadpLol + {z - zo)dpoJa) dx p , (4.20) 

and consequently, from the definitions of the Frank and Burgers vectors, 

lim n!. = < {earBl - {z - zo)Q,l, - xoMl}daSo, >, (4-21) 

£-»0 + 

where is the 2D Dirac measure located at and tpz{xa) := tp{xa,z), while symbol 
<C •, ■ 3> denotes the 2D distribution by test-function product. 
• When k — K, definition (32] shows that 



dpbl := + e^T(xj — X0'y)di3Uj*, 
from which, after multiplication by e^a, it results that: 

XadrlJ^K ~ —(.Kadrhl + e^a£^z + XOadrl^% + {Xk — X0K)dTOJa- 
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Then, by lemma IT2l with a permutation of n and a, (|4.19|) also writes as 
n'„ = daf>{x) I {-e^adpbl + eKaS*ftz + xoadpuil, - xondpuja) dxp + 0(1). 

On the other hand, from equation (|4.19p and lemma|42](i.e. from strain compatibility) 
it follows that: 



n'„ = dc,ip{x) [-ei^pSffzdXa + eapSfizdx^) + o{l) 

= dcip{x) [ e^^£lzdx0+o{l). (4.22) 



By summing this latter expression of 11^ with (I4.22|) , from the definitions of the Frank 
and Burgers vector it follows that 

n', = ia<,<^(i)e<.« {Bl - e^ipn^-^xop) + o(l). (4.23) 

Hence, in the limit e — >■ 0^ (|4.23|l writes as 



lim n'^ =< \ -e^aB* - -eKae^0^*xoi3 } daSo,<pz > ■ (4.24) 



;-o+— " ^2 2 
▲ Therefore, the result is proved on Q,^, since 

lim nfc(2,(^,e) = lim n'fc(2, </?,£)+< nfc(5o,</3z > . (4.25) 

e->0+ £^0 + 

As suggested by equation (|4.13p . to obtain the result for the entire domain Cl it suffices 
to integrate equations (|4.20|) and 1)4. 23p and the expression QXip{x) over L, in order to 
replace So by the line measure 5l in (|4.21l) . (|4.24|) and (|4.25p . By (|4.13|l the proof is 
then achieved. □ 

Theorem 4.2 (main result for a set of isolated defect lines) Let in the2D case 
L'*', i G I C N stand for a set of isolated parallel dislocations and/or disclinations 
passing by {x''p\ z) and fiz'*' , _B^*'^ and Sj^(i) denote the associated Frank and Burgers 
vectors, and the concentrated ID Hausdorff measure density on L'*'. Then under as- 
sumptions \41\ o,nd \4.S\ in the vicinity of each defect line, incompatibility develops as 
the distribution 

Vk = Skzvl + SkK-nl, (4.26) 

with 

r,: = + (^7*'^ + - xofi)n:^'^) d^S^u) ) , (4.27) 

Vl = ie.,^B:«9,5^(„. (4.28) 

Proof, the proof is straightforward from theorem 14.11 An alternative formulation is 
provided by ^^-^S^. □ 
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4.3 Applications of the main result 

Throughout this section, {x,y,z) denotes a generic point of Ql and ah tensors are 

written in matrix form in the Cartesian base {ex,ey,ez). 

• Screw disclocation. Since B*^ = ill = 0, P-llf) and (|4.12|l yield 



2 



dySL 





This result is easily verified with use of equation p.Sp . One needs to compute < 
77^,</9 / ekpn£ap£^„dpdc,'pdV , that is to calculate the integral of 



4-K 







By integration by parts, using Gauss-Green's theorem on fl, and recalling that test- 
functions have compact supports and that dm logr = — ^, this integral becomes 



./n 



dy^{d.^+dy'-^) 





dV - 



47r 



-dyipdm^ogr 




dV. 



Hence, from the relation A (logr) = 2n5L, the first statement is verified. 

• Edge dislocation. Whereas dmijJk identically vanishes on it is easily seen that 
||TTT|| and (jn^ with B* = SI* = yield 



[vl] = b; 








We must compute < r]k,ip >— / ekpn^aisS/sndpda'fidV . For n 7^ 3, the strain com- 

Jn 

ponents do not identically vanish and, for k = 1 and k — 2, we must have p = 3 
and hence the only non- vanishing component of the expression ec,0£0„dpdaf are 
£yxdzdyip — Eyydzdxifi soid Sxydzdxifi — Sxxdzdyip. By integration by parts, recall- 
ing that the strain does not depend on z, the related integrals vanish. For fc = 3, the 
integrand is 

Cpnzf^apSlndpda^ = {dySxx - dxSxy) dylfi + {OySxy - dxSyy) dxlfi. 

By inserting the expression of the strain tensor into the right-hand side of this equation, 

integration by parts provides the expression / -dx'pA{\ogr)dV, achieving the 

in Stt 

second verification. 



Wedge discUnation. Incompatibility reads 

[vi] = ^: 






5l 
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rjl,ip>= / ekp„€a0£^ndpda<fidV. For k 
Jn 



We must calculate < rjl, (p >— / €kpn^a$£gndpda(pdV . For k — I and k = 2, we must 

In 

have n 7^ 3 and p = 3, but then the integrand vanishes. For A; = 3, we compute 
ep„ei„^£mndpdif = ((9A(log— )H — (pA{log—) 

achieving the third verification. 

5 Conclusive remarks 

In this paper a general theory revisiting the work of Kroner (1980) has been developed 
to model line defects in single crystals at the mesoscopic scale. A rigorous definition 
of the dislocation and disclination density tensors as concentrated effects on the defect 
lines has been provided in the framework of the distribution theory. The main difficulty 
resulting from the multivaluedness of the displacement and rotation vector fields in 
defective crystals has been addressed by defining the single-valued Burgers and Frank 
tensors from the distributional strain gradient. Whereas outside the defective lines 
both tensors are regular functions directly related to the displacement and rotation 
gradients, in addition they exhibit concentrated properties within the defect lines 
which may be linked to the displacement and rotation jumps around these lines. 

Moreover, defining the incompatibility tensor as the distributional curl of the Frank 
tensor, the principal result of our work has been to express in the two-dimensional 
case incompatibility as a function of the dislocation and disclination density tensors 
and their distributional gradients, and to demonstrate this relationship under precise 
assumptions on the regularity of the strain tensor in the vicinity of the assumed isolated 
defect lines. In a subsequent paper (Van Goethem & Dupret, 2009fo), our theory is 
extended to the case of a countable number of defect lines under specific hypotheses 
based on the geometric measure theory. 

In general our work is devoted to provide a rigorous distributional definition and a 
new understanding of the different mathematical objects (dislocation and disclination 
densities, contortion, incompatibility. Burgers and Frank tensors, elastic strain, etc.) 
that can be added at the mesoscopic scale in order to well-define the associated ho- 
mogenized objects at the macroscopic scale. Further work will deal with the general 
three-dimensional dynamic theory. 
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Abstract 

This paper develops a geometrical model of dislocations and discli- 
nations in single crystals at the mesoscopic scale. In the continuation 
of previous work the distribution theory is used to represent concen- 
trated effects in the defect lines which in turn form the branching lines of 
the multiple-valued elastic displacement and rotation fields. Fundamen- 
tal identities relating the incompatibility tensor to the dislocation and 
disclination densities are proved in the case of countably many parallel 
defect lines, under global 2D strain assumptions relying on the geometric 
measure theory. Our theory provides the appropriate objective internal 
variables and the required mathematical framework for a rigorous homog- 
enization from mesoscopic to macroscopic scale. 

1 Introduction 

1.1 Preliminaries and principal hypotheses 

The present paper provides a mathematical theory of the geometry of crystal disloca- 
tions and disclinations in continuation of the work of Van Goethem & Dupret (2009a) 
where the general context of this research is detailed and which will be referred to as 
Part A in the sequel. In summary, the objective of these investigations is to develop a 
rigorous mathematical framework for the treatment of line defects in single crystals at 
the mesoscopic scale. As explained in Part A, concentrated effects in the defect lines 
and their neighbourhood have to be considered at this scale and this requires to make 



1 



use of the distribution theory (Schwartz 1957) to handle the related fields (dislocation 
and disclination densities, contortion, incompatibility, etc.) and their relationships. 
Moreover, in view of the incompatibility of the elastic strain tensor in the presence of 
line defects, the associated rotation and displacement are multiple-valued fields whose 
branching lines are precisely the defect lines. The combined treatment of distributions 
and multivalued functions was addressed in Part A, where our theory was applied to 
the case of a set of isolated parallel, moving or not, line defects under the hypothesis 
of a 2D elastic strain field. 

In this second paper, the case of countably many parallel defect lines is investigated. 
Therefore, instead of analysing the regularity of the elastic strain near an assumed 
isolated defect line, a more general abstract approach is selected with a view to defining 
the appropriate functional space to validate the main theorem relating the strain 
incompatibility to the defect densities. 

Let us remark that our mesoscopic setting will be able to treat fine and complex 
dislocation structures since accumulation lines or points in the defective set will be 
allowed (such as typically the structures appearing in the work of Cantor (1915) on 
transfinite numbers, see figure [T]). This feature represents a key ingredient of our 
theory since a tending to infinity number of defect lines unavoidably appears in the 
homogenization process from meso- to macro-scale. Moreover, when the defect lines 
exhibit a clustered mesoscopic structure, even if their actual number will always remain 
finite across any bounded area, it is much more convenient mathematically to consider 
a model where this structure may be infinitely refined. 



Figure 1: uP dislocation structure, where u is the first ordinal transfinite number 
(Cantor 1915). 

Before any further development, precise definitions and assumptions are required. 
In general, the functional spaces used consist of distributions. Radon measures (Am- 
brosio et al. 2000, Mattila 1995), functions, etc., which can be considered as continuous 
functionals over a set of test-functions whose regularity determines the functional space 
properties. However particular care has to be given to avoiding undesirable boundary 
effects. 

Definition 1.1 (functional spaces used) The crystal domain is an open connected 
set on which some mathematical elements ( distributions, Radon measures, locally 
summable functions, etc.) are defined as Imear functionals over an associated set of 
test-functions whose support is a compact subset of il. Henceforth, the qualification 
"on Q" for these elements will mean in addition that extensions of these elements 
exist as functionals over all the test-functions having the desired properties and whose 
support is a compact subset ofW^. 

In other words, a distribution. Radon measure, locally summable function on Q 
will always be constrained to also be the restriction to S7 of a mathematical element 
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of the same type defined on the whole R'^. A simple ID example can illustate this 
constraint. The function 

/ : R,| — ^ R, XI — y f{x) = 1/x {x > 0) (1.1) 

is locally summable on Rj in the classical sense but cannot be extended as a locally 
summable function over R, this resulting from its behaviour for x ^ 0^ . According to 
definition ll.il in the present paper this function will not be called locally summable 
on . Indeed, in a context where locally summable functions are needed and where 
the physical domain of interest is R,"}^, no shift to the right of / as defined by (|1.1|) 
can provide a locally summable function on Rq" (with g{x) = f{x — a), x > a > 0) 
whatever the definition of g is for < a; < a. So, if to be acceptable the locally 
summable functions on Rj are requested to exhibit the same properties near the 
origin as in the vicinity of their interior points (to avoid peculiar boundary effects for 
x — ^ 0"*"), freely translating these functions to move the origin onto an interior point 
should be allowed, and hence / cannot be accepted as locally summable on R|J" if g 
cannot. This justifies definition ll.il 

Similarly, the same function / defined by (|1.1[) can be considered as the density 
of a Radon measure on Rj in the classical sense (as acting against continuous test- 
functions whose support is compact and contained in Rq"). However, again no shift to 
the right of / can provide a Radon measure on RJ since no extension to R of / as a 
Radon measure in the classical sense exists. So, in our work / will not be considered 
as a Radon measure on Rq" according to definition ll.il 

Considering now / as a distribution on Rj (as acting against C°° test-functions 
whose support is a compact subset of Rq"), / can be prolonged as a distribution on R 
by defining for example the prolonged / as the pseudo- function Pf.{l/x), which is the 
distributional derivative of the function F{x) := log |a;| and acts again a test-function 
by taking the Hadamard finite part of the resulting diverging integral. So, here / can 
be considered as a distribution on RJ in the sense of definition ll.il 

It should be noted that the spaces generated from definition 11.11 are non-closed 
subspaces of classical spaces (the usual distributions. Radon measures, ... on Q) and 
hence cannot share all their properties (completeness, etc.). 

In the crystal domain fi, the meso-scale physics will then be represented by a 
nowhere dense set of defect lines which in 2D are parallel to each other. 

Definition 1.2 (2D mesoscopic defect lines) At the meso-scale, a 2D set of dis- 
locations and/or disclinations £ C is a closed set of f2 (this meaning the inter- 
section with ^ of a closed set of R'^ j formed by a countable union of parallel lines 
L^^\i £ I C N, whose adherence is itself a countable union of lines and where the 
linear elastic strain is singular. In the sequel, these lines will be assumed as parallel to 
the z-axis. 

The present mesoscopic theory will be completely developed from the sole linear 
strain - which itself could be defined from the stress field (although the stress-strain 
relationship is not used in the sequel) and therefore is an objective internal field. 

Assumption 1.1 {2D mesoscopic elastic strain) The linear strain E^^i is a given 
symmetric L\^^{Q,) tensor (in the sense of deHnition \l. 1\) prolonged by on the dislo- 
cation set C, and such that d,£^n = 0. Moreover, is assumed as compatible on 
0,c 0,\ C in the sense that the incompatibility tensor defined by 

INCOMPATIBILITY: rjli ■- etp^eig^dpOgS*,^, (1.2) 
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where derivation is intended in the distribution sense, vanishes everywhere on Qc- 
Equivalently, in 2D there are real numbers K, aa and b such that the following equalities 
hold on ^Ic: 

ea-^eissdadpE^s = 0, 

eapdaS}^ = K, (1.3) 

£zz — OaXa ~\~ b. 

In Part A, key tensor fields were derived from the mesoscopic elastic strain field 
as order-1 distributions (acting on test-functions with compact support). 

Definition 1.3 (Frank and Burgers tensors) 

FRANK TENSOR: dm^^l — ekpqdpS*^, (1.4) 

BURGERS TENSOR: dibl := + ekpq{Xp - Xop)diUJ*. (1.5) 

Line integration of the Frank and Burgers tensors in Qc (i.e., outside the defect set) 
provides the index- 1 multivalued rotation and Burgers vector fields ljI and bl (with 
"index-1" meaning that their first derivatives in Qc, denoted by S^^'cjJ and ^j^'fo^, 
are single- valued) . These properties are summarized in the following theorem, whose 
proof is classical. 

Tlieorem 1.1 (multiple- valued displacement field) From a symmetric smooth 
linear strain Eij on Qc and a point xo where displacement and rotation are given, 
a multivalued displacement field u* of index 2 (whose second derivatives are single- 
valued) can be constructed on Qc such that the symmetric part of the distortion d^^'^u* 
is the single-valued strain tensor f * while its skew-symmetric part is the multivalued 
rotation tensor uj*j ;= — eijfecj^. Moreover, mside Qc the gradient Sj"' of the rotation 
and Burgers fields uj^. and b'^ — — ekim^i {xm — xom) coincides with the Frank and 
Burgers tensors. 

From this result, the Frank and Burgers vectors can be defined as invariants of any 
isolated defect line i'*' of £. 

Definition 1.4 (Franlc and Burgers vectors) The Frank vector of the isolated de- 
fect line L'*' is the invariant 

Qf'>:^Hf\ (1.6) 

while its Burgers vector is the invariant 

Bf^ ~ [blf^ = [ul]^'Hx) - ekir.Q;^'Hx^ - xo^), (1.7) 

with [ajfe]''-*, [fefe]''-* and [u%]^''^ denoting the jumps ofujl,bl and u^ around L^^\ 

The case of non-isolated defect-lines represents a major issue of this work and will 
be resolved at a later stage. Besides their relationship with the multivalued rotation. 
Burgers and displacement fields, the Frank and Burgers tensors can be directly related 
to the strain incompatibility by use of p.2|l . (|1.4p and p.5|l . 

Theorem 1.2 The distributional curls of the Frank and Burgers tensors are 

tiijdidjujl = 77*fc, (1.8) 
ea^didjbl = ekpq{xp - xop)riiq, (1-9) 
with rj*^. the incompatibility tensor. 
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From this theorem it results that single- valued rotation and Burgers fields and 
can be integrated on Q if, and only if the incompatibility tensor vanishes. 

The dislocation and disclination densities are the basic physical tools used to model 
defect density at the meso-scale. In 2D (here meaning that dz£mn ~ 0) the defect 
lines are parallel to the z-axis and hence only some components of the defect densities 
do not vanish. 

Definition 1.5 (2D defect densities Q) 

DISCLINATION DENSITY: 0* := ^ , (1-10) 

DISLOCATION DENSITY: Al := ^ Bl'-"^S^(i) (fc = 1 ■ ■ ■ 3). (l-H) 

In general, the complete defect density tensors are denoted by Q*j and A*j with the 
indices i and j referring to the local defect line direction and Frank or Burgers vector, 
respectively. So in 2D, 

eL = olk^o, eL = e: (1.12) 
a:^ = 0, a:, = A^ (1.13) 

with fc = 1 ■ • • 3 and k = 1, 2. The vanishing of OJ„ was shown in Part A. 

An additional important defect density tensor called the contortion was introduced 
by Nye (1953) from the work of Kondo (1952). 

Definition 1.6 (2D mesoscopic contortion) 

CONTORTION: K*j := S^^a* - ^a^Sij (i, j = 1 • • • 3), (1-14) 

where 

al := Al - Ska^aiiQlixis - xois)- (1-15) 

Here, al is an auxiliary defect density measure associated with a general 3-D tensor 
a*j = A*j + ejimQ*ii{xm - xom) such that in 2D 

Q*fe = 0, Q*fe=afc, (1.16) 

with k = 1 ■ • • 3 and k = 1, 2. The general 3D contortion is k*j — a*j — ^alj^„^Sij, in 
such a way that in 2D, k*^- vanishes if a 7^ j. 

I. 2 Objective of this work 

The principal objective of this work is to identify key distributional fields at the meso- 
scopic scale and to demonstrate their relationship in a rigorous functional analysis 
context, with a further view to providing the required framework for the homogeniza- 
tion of these fields and their relations to the macroscopic scale. 

In 321 the main theorem of Part A (expressing the elastic strain incompatibility 
in terms of the defect densities and their gradients) will be extended to the case of a 
countable ensemble of parallel defect lines. To this end, besides the strain assumption 

II. ll an additional assumption is made on the Frank tensor p.4|l . 
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Assumption 1.2 (mesoscopic strain assumption) The 2D strain belongs to 
L}^^{i^) (in the sense of definition \l.l\) and is compatible on flc- Moreover, the {m,z) 
components of the Frank tensor dm^* = eapdaS^rn (1 3) are Radon measures 

on O,, whose singular parts with respect to the Lebsgue measure are purely concentrated 
on £ while their absolutely continuous parts have a 2D curl which itself is a Radon 
measure on Q. 

Remark 1.1 No similar assumption could be made on the complete Frank tensor 
without contradicting the edge and screw dislocation examples of Part A. Moreover 
the absolutely continuous part of dmij-'t cannot be required to be of bounded variation 
without contradicting the wedge disclination example of Part A. On the other hand, 
it will be seen that the sharp assumption \1.S\ is required to establish our theory in the 
general case of countahly many dislocations. 

Among several equivalent formulations, our main theorem then takes the following 
form. 

Main theorem: 

incompatibility decomposition for a countable set of 2D defect lines. 

Under assumptions \1.1\ and \l.Sl incompatibility as defined by SI. 2}) is the vectorial 
first order distribution 

V*k = Vki = 9*fc + uijdinlj. (1-17) 



Let us also introduce here the main intermediate proposition needed for the proof 
of the above representation theorem. Under assumptions 11.11 and 11.21 the following 
decomposition theorem will be proved in the 2D case. 

Theorem: 2D strain decomposition. 

Let the 2D strain tensor be a compatible L\^^ tensor on Q.c . Then the following 
decomposition holds: 

where e^,i is everywhere compatible, whereas is a sum, 

E* = (1 19) 

where each Emn {i € T) is analytically known, compatible and smooth on Q,\ L^*', 
while E^„ is singular on L and compatible and smooth on ^Ic- 

From our main theorem, i|3]will then be devoted to introducing new mesoscopic 
distributional fields, namely the completed Frank and Burgers tensors, which will rep- 
resent the appropriate objective internal variables after homogenization to the macro- 
scale, and to reformulate the main theorem in their terms. Conclusions will be drawn 
in g] 
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2 Distributional analysis of incompatibility for 
a countable set of parallel dislocations 

To capture the macro-scale physics, homogenization must be performed on a set of 
dislocation lines whose number tends to infinity in order to define diffuse defect density 
tensors. Therefore, assumption 11.21 was introduced in a functional formulation that 
can be extended in some way from a set of defect lines (at the mesoscopic level) to a 
diffuse defect density (at the macroscopic level). 

The extension of our theory to a countable number of defect lines poses several 
technical problems. A first difficulty arises from the different kinds of convergence that 
could be required. Typically, considering a series of Dirac masses on /q = £r]{z = zq}, 
then its convergence as a measure implies that the sum of the weights must converge 
absolutely, but this is no longer the case if a (coarser) distributional convergence is 
required. A second example is provided by those distributions that are the gradient of 
a summable function. If these distributions are concentrated on isolated points, they 
must be the sum of Dirac masses, whereas this property might fail on a countable set. 
More generally, it is known (Schwartz 1957) that a concentrated first-order distribution 
on isolated points is a sum of weighted Dirac masses and Dirac mass derivatives, while 
a concentrated measure on a countable set is a sum of weighted Dirac masses. However, 
it is false to claim that a concentrated first-order distribution on a countable set is a 
sum of Dirac masses and Dirac mass derivatives, as ID counter-examples can show. 
In general, a more complex mathematics governs the accumulation points of Iq, and 
appropriate tools are required to extend the representation theorems of Part A to a 
countable set of defects. 

2.1 General strain decomposition property 

In general any vector field can be decomposed into a solenoidal and an irrotational part, 
and this property can be easily extended to distributional fields. In this paper, the 
similar decomposition of any symmetric tensor field into a compatible and a solenoidal 
part will be used to extend the main theorem of Part A from isolated to countable 
dislocations 0. In what follows, we will first give a proof of the decomposition existence 
in the general distributional case and then investigate its regularity in the 2D case. 
The main theorem will be extended in a further section. 

Theorem 2.1 (standard decomposition of a symmetric tensor) Any symmet- 
ric distribution tensor E (or Eij) can be decomposed into a compatible and 
a solenoidal symmetric part: 

E^E^ + E", (2.1) 

with 

V X E" xV = (compatible E"), (2.2) 

and 

V ■ E" =0 (solenoidal E"). (2.3) 

^Kroner (1980) first observed tliat this decomposition provides a link between the disloca- 
tion density in a medium and the associated strain tensor incompatibility. 
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Proof. A Any tensor E" defined by the relation 

= V X F X V (2.4) 

is symmetric and solenoidal if F is a symmetric tensor distribution. Then the reminder 
E'^ — E — is compatible provided, after some calculations, F satisfies the relation 

AAF,, + d.djdkdiFu ~ A (djdkEk) - A (dAF.k) = e,kiejmndkd^Ei„, (2.5) 

with A the Laplacian operator (A — didi). If in addition the gauge condition 

V • F = (2.6) 

is imposed, then 1)2. 4|) reduces to the elliptic equation 

AAF = V X F X V. (2.7) 

▲ Therefore, to find the searched decomposition (|2.1|) . (|2.2|) . (|2.3|) . it is sufficient 
to solve 1)2. 7|) for F with the gauge condition (|2.6|l . If E is sufficiently regular, F 
will simply be found by solving (|2.7p with, among others, the 6 boundary conditions 
\7 ■ F = and d {V ■ F) /dn = 0. As a matter of fact, a solution exists because 
the operator AA is elliptic, and this solution is divergence-free because taking the 
divergence of 1)2. 7|l provides the relation AA (V ■ F) = which, together with the 
boundary conditions implies that V • F itself vanishes. 

If E is not sufficiently regular, E can be approximated as a distribution by a family 
of C°° functions F,(e > 0) with F, ^ F for e 0+ (Schwartz 1957). The family of 
equations obtained by replacing F by Fe in (|2.6[) . (|2.7p provides a family of solutions 
Fe which tends to a suitable F when e — >■ 0'*'. □ 



2.2 First representation theorem of a 2D incompatible 
strain 

The previous section has shown that a distributional decomposition of the symmetric 
strain £* G L}a^{Q,) into compatible and solenoidal distributional parts £*'^ and £*" 
always exists, with the right-hand side of (|2.7p showing to be the incompatibility 
tensor. However, more regular solutions exist in the 2D case. Before proving them, 
the following result will be needed. 

Lemma 2.1 Let 5''' stand for the Dirac measure at x''' G lo and ^^C'''^''' be a 

Radon measure on S7,(, = Q, \ {z = zq} in the sense of definition \l.l\ Then the sum of 
the weights C'*' is locally absolutely convergent, this meaning its absolute convergence 
on any bounded subset £ X'} oflo. 

Proof. Since C'''(S''' is a Radon measure, then C'''^'*' is a finite Radon 

measure and the sum can be indifferently carried out on every permutation of T' . 
Hence, taking a test-function wich equals 1 on io, the sum of the weights converges for 
every permutation of I' and is absolutely convergent. □ 



8 



Remark 2.1 If^C^''>6^'^ were assumed to be a general distribution instead of a 

Radon measure, no such statement on the absolute convergence of the sum of the 
weights could be proved as the following simple ID counter-example shows: selecting 
x'-^^ = 1/i {i G 2: = No) andC'-^'' = (-1)'+^ + 1/(1 + i)) provides a distributionally 
convergent series C''';?'''', since it is the derivative of the L\^^ converging series 

iGI 

1 - ^C*''(l - //*'^) with H^^ = H{x - x''^) and H the step function, whereas the 
sum |C'''| does not converge. 



Notations 2.1 Henceforth {a;''', i G X} will denote the set of points defining Iq. 

Theorem 2.2 (regularity of the strain decomposition) Let the strain and the 
Frank tensor satisfy assumptions and M.'A and the dislocation set be defined ac- 
cording to definition \l.S\ Then the following decomposition holds: 

where £ is compatible, whereas £ Ll^^{Q,) is solenoidal. 

Proof. Consider any 2D cut Q^g of SI and assume first that fl^g is bounded (extension 
to unbounded sets is direct). Since the strain is independent of z, in 2D it suffices to 
solve (|2.7|l and (|2.6|l for F on SI^q with Eij = £*j. This will be achieved by solving 
an associated problem on Q^o by means of complex (but not necessarily analytic) 
functions of two real variables. To this end, (|2.7|) is first expressed in block matrix 
notation: 



" AAFcp 


AAFcz ' 








_ AAFz^ 


AAF 




via 


viz 



d^£* 

-dxdy£zz 


-dxdy£l, 
dl£lz 


-dx {dx£yz - 


dy ^■xz) 
- dy£xz) 


dy [dx£yz 
~dy£xz) 


{^x£yz 

~^y^xz) 


dx i^^X^yy 
dy {px^xy ~ 


dy£xy) 
dy£xx^ 



(2.9) 



and every block equation is separately solved. 

▲ First block (rj^p). By the compatibility condition (|1.3p outside Iq, it results that £*^ 
is linear on the open and connected defect-free region Qzg \ lo and can be prolonged 
by a linear function on Qzg- Since Fai3 = is a solution of AAF^^ = on Qzq, by 
112.911 an admissible £tz is 



£: 



0. 



(2.10) 

A Second block (rjaz)- As dzUJ^ = dx£yz ~ dy£xz, it is convenient to seek a solution of 



AA {Fxz + iFyz) = -i [dx + idy) dzOJ^. 
According to assumption 11.21 dzio* decomposes as follows: 

iex 



(2.11) 



(2.12) 
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where the absolutely continuous part of the now finite Radon measure dzLU* shows to be 
the constant K appearing in the compatibility condition (|1.3p . while its singular part 
is purely concentrated. Moreover, the sum of the weights c'*' is absolutely convergent 
by lemma [2Tl1 when Qz„ is bounded. So (|2.1H) develops as 

AA (F,, + iFy,) = -i (a. + idy) (2.13) 

Then, since A A rewrites as {dx + idy)^ {dx — idy)'^ , it suffices to solve 

{dx - idy) {Fx. + iFy,) = -iJ^, (2.14) 

with T a solution of 

= Yc^'^S^'\ (2.15) 
To solve this system, observe that (|2.14|) develops as 



QxFyz OyFxz , 

dxFxz QyFyz — 0, 

where an acceptable W^'^ {Q.zo) field T satisfying (|2.15|) is 



(2.16) 



using the notations x = (2^,2/), '0 = {x'*' = (2;'*', 2/'*'), i G 1} and r'*^ := |x - x'^'l- 
Hence, by p.4|l an admissible solenoidal strain belonging to L^{Q.zo) is 

£Tz = dy{dxFyz~dyFxz) = ~dyT, (2.18) 

C = 'dx{dxFyz-dyFx,)^dxT. (2.19) 
▲ Third block {rj*^). Recall first that 9^ daS*fp and that 

V*zz = 9a [daSlf) - d0£*0) = ea/jdad/lUll, 

in such a way that the problem is to solve 

AAF = dxdyLo*, — dydxuj* = ^{i{dx + idy){dxUJz — idyujl)} (2.20) 
where, according to assumption 11.21 daCd* develops as 

a,c.: = ^c«5<''+/. (2.21) 

with absolutely convergent sums of the weights c^' by lemma ETT] when flzo is bounded 
and for some functions whose curl is here a finite Radon measure, which must be 
concentrated on lo to ensure the compatibility of outside the defect set. So, (|2.20|l 
rewrites as 



AAF 



{dx + idy) [ ^ (cW - icW) 5'') J I + dxfy - dyf. 
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Now, in view of the properties of fa resulting from assumption 1 1.2 1 the last terms 
write as 

=^C«5«, (2.22) 

where the sum of the weights is absolutely convergent. Eventually, using the distribu- 
tional identity 



with the notation 5a;''' = x — x^'\ Sy^'^ ~ y — y^''\ (|2.20|) can be written in the form 
AAF = K + idy) (j2 (c« - <5« - Y: ^ {y'-' + • 



(2.23) 



A particular solution of (|2.23l) is provided by solving 

HO 



(9. + idy){d. - idyf{F + iH) = i(Yl (c-^ - ^4'') 

iei / 
with H an additional unknown. This latter equation is equivalent to the system 

(d^ - idy){F + iH) = g on (2.25) 
AQ = i(dxid* — idyuj*) on Qzg, (2.26) 

which can be easily solved. In a first step, a particular solution of (|2.26|l is given by 
g = ^1+52, 



^2 = (2.28) 



with both Qi and C/2 belonging to W^'^{Q,g). In a second step, (|2.25|l is simply 
rewritten as 

d^F + dyH = ^{g} and d^H - dyF ^ ^{0}, (2.29) 

whose solution F = Fi + F2 and H = Hi + H2 is given by 

Fi = dxi>i + dyifii and Hi = dyiii — dx^fi, (2.30) 

F2 = dxi>2 + dy(p2 and H2 = dyii2 — dx^2, (2-31) 

for some gauge fields Tl}i,(pi,^2,'P2 satisfying the equations 

At/)! = and A(/pi = (2.32) 

A1P2 = K{g2} and Aifi2 = -3{e2}. (2.33) 
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Particular solutions of (|2.32|) belonging to W^'^ (Q^q) are 



i)2 



i6I 

in such a way that 



log rt'> - 1 
8^ 



and ifi = — ( 



(i) 



Stt 



i6I 



belongs to W'''^{Qzg), thereby defining the solenoidal strain f 



(2.34) 



(2.35) 



which belongs to L^{Qzg) r\C°°{flzo \ lo)- Similarly, particular solutions of p.33|) 
belonging to W^'^{Qzo) are given by 



1p2 



'P2 = 



327V 



and hence 



^ ^ 327r ' 



F2(a;,y) = ^^r«^ (21ogr(') - 1)) 



also belongs to W^'^{^lzo), defining the solenoidal strain f*^''^: 



dlF2 ~d^dyF2 



(2.36) 
(2.37) 

(2.38) 
(2.39) 



~d^dyF2 d'iF2 
which belongs to L^i^zg) n C°°{nzo \ lo)- 

▲ Summary. The solenoidal part of the strain is the tensor £^„, 

^mn — £ zz^^^z^nz ~\~ ^az i^^mz^na ~\~ ^ma^nz) ~t~ (*^a/3 ^ajS ) ^^ct^n^ ^ (2.40) 

where f^l, C^'^ and f*^'^ are given by (f2T^ . (pTT]) . pH]) . (p^Sj) and ((2391), 

respectively, and all belong to L^{flzo), observing that the weig hts c('',cL'' and C''^ 
defining the intermediate expressions T,F\ and F2 in (|2.15l) . (|2.34|) and (|2.38|l are 
associated with absolutely convergent series: 



^ic("'|<oo, ^||c<:'|i<oo and ^|C<*'|<oo. 



(2.41) 



Therefore belongs to L^{Q,zq), is solenoidal and satisfies (|2.8p . Extension of the 
proof to unbounded domains Q.z„ is immediate. □ 
The coefficients c^'' , ci* ' and C^'' will show to be the Burgers and Frank vectors of 
screw and edge dislocations and wedge disclinations, respectively. 
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Remark 2.2 The hypothesis provided by assumvtion \l.lA that dmUJ* has an absolutely 
continuous part whose curl is a Radon measure is a request to make the proof in the case 
of a countable set of line defects. Indeed, when the 2D defect set Iq has accumulation 
points in Q^o , a complex distributional behaviour can take place near these points which 
forbids geting the proof if a sufficiently strong hypothesis is not introduced to account 
for a possibly countable number of disclinations on the sole basis of the strain field 
properties. More tractable hypotheses on daUj* itself (and not its curl) are currently 
under investigation. 

As a ID example to illustrate the above difficulty, the function 

F= J2 C-'^H" - H'-'^) 

with H^'^ = H{x - x^''>), x^''> = 1/i, H° = H{x) and H the step function, may 
correspond to an L}^^ converging series even if the sum of the weig hts C<'' diverges. To 
show this, it suffices to select appropriate C'*' such that the partial sums defining F are 
all enclosed between the Lj^^ functions G{x) and —G{x), with G{x) = log ((1 + a:) /x) 
for X > and G{x) = for a: < 0. The Lebesgue dominated convergence theorem 
then shows that F £ L}^^, in such a way that the distributional derivative of F, 
which cannot be the diverging series — C'*'^'''' , exhibits a special behaviour near 

iei 

the origin to recover convergence. Similar effects take place in 2D and appropriate 
assumptions are then necessary to obtain (|2.22|) . 



2.3 Second representation theorem of a 2D incompatibile 
strain 

This section provides a further decomposition of the strain, since the solenoidal part 
is itself decomposed into an everywhere compatible part and another smooth part 
outside from the defect set £,. 

Theorem 2.3 (analysis of the singular part of the strain decomposition) Let 

the strain and the Frank tensor satisfy assumptions and \1.S\ and the dislocation 
set be defined according to deHnition M.Sl Then the solenoidal component of the strain 
satisfies the following decomposition: 

^mn ^mn -^mn ; (2.42) 

where S^n G iLc(^zo) ^■^ compatible on and where 

E*m^ = e ^'o- (") ' (2-43) 

with Emnii £ T) smooth and compatible onQ, \ L^''. Moreover, the Frank tensor part 
tkpndpE^n is smooth on Q,c. 

Proof. As in the previous proof, since the strain is independent of 2, it suffices to work 
on the 2D domain Q.^^ which is again assumed to be bounded without loss of generality. 



▲ £a~ components. Following decomposition (|2.8p . write 

£Tz ~ dy (dxFyz — dyFxz) — E*^ + S'^z, 
£yZ (^dxFyz dy Fxz) Eyz ^ ^ y 
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with vanishing E'^f and £y 



^xz yz 



Hence, according to (|2.18p and (|2.19p . E'^^ and E*^ are sums of screw dislocations: 



(2.44) 



which from ((2^ belong to L^{^zo) n C°°{Q.zo \ The relations (fOi)) also show 
the smoothness of the Frank tensor parts eapdaE*^ and e/^jd^E^-, outside Ig. 



A £11 and components. By (|2.10p the expression 

E* := Sit = 



(2.45) 



exhibits the form (|2.43p and provides Frank tensor parts ea-yd^E*^ which identically 
vanish on Qzq- On the other hand, it can be checked that on Q^o \ io'- 



E 



5:rW(l-2^) -V')(l-2^) 
-5y«(l-2^) 5^«(l + 2^) 

5y'''(l + 2^) -5^«(l-2^) 

-5:rW(l-2^) 5y«(l-2^) 



(2.46) 



/*c,l 



Define then in the decomposition (|2.42p . (|2.43|l of f*^'^ the components E*^p and f'*'^ 
as follows: 



(2.47) 



^E 
E 



4j^j.(i)4 



4j^j.(i)4 



5y(>)(52;(.)2 _ 52/(>)2) fe(«)(52;(>)2 + 35y(i)2) ■ 



(2.48) 



Calculations show that S'^p^ is the difference between (|2.46p and (|2.47p . is compatible, 
and also belongs to L^{^lzo) by (|2.4ip . Moreover E^^^^ is a sum of edge dislocations, 
whose curl is vanishing and which are smooth on fi^^ \ {(i" , y «)}. 



▲ f components. Define 



1-1*2 r'*s,2 



with vanishing f Calculations show that 



Q/3 



47r 



(2.49) 



(2.50) 



in such a way that is a sum of wedge disclinations, which are smooth and of 
vanishing curl on Vlzo \ {(£'*', y*'')} (cf. Part A). 
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▲ In summary, the following solenoidal strain decomposition has been proved: 



where 

^mn • — ^zz ^mz^nz ~t~ ^az {^^ma^nz ~\~ ^mz^na) ~\~ (^^^ ^afS 

is compatible on fi. Moreover, writes as 



n0 



(2.51) 



(2.52) 



(2.53) 



with E*^,E*p^,E*} and E^} defined by 12351), (p:T7|) and (^301) . and hence 

Emn together with the related Frank tensor ekpndpE^„ is smooth on Qc thereby 
terminating the proof. Extension to unbounded sets ^Izq is straightforward. □ 

2.4 Applications of the strain decomposition 

In this i]2E.4l X' refers to any bounded subset {i'*', i G I'} of Iq. 

Set of parallel screw disclocations. Part A and (|2.44p directly provide the equality 



_B*''' and a vanishing strain compatible part. Moreover according to (|2.4ip the 



following condition holds: 



^|b:«|<oo. 



iei' 



Set of parallel edge disclocations. From Part A and (|2.46|) - (|2.48l) . it turns out that 



By and Cxz = , with according to (|2.41|) the following bounds: 
^|B:'"'i<oo and ^|B;<''|<oo. 



Set of parallel wedge disclinations. The expression of [£ij] given in Part A : 



4-K 



1 + log r 
1 + log r 




n*(i + t^*) 

Stt 



cos 261 
sin 26 




sin 261 
- cos 261 




which shows to be the sum of 



log r + sin 6 — sin 6 cos 9 
— sin 6 cos 9 log r + cos^ 9 




and a compatible part (since the associated Frank tensor part vanishes). Therefore, 
according to (g^DJ, C''> = Q*z^'^ with by f^Ul} the bounds 



< CO. 
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2.5 Mesoscopic defect densities in 2D incompatible elas- 
ticity 

The following theorem expresses the 2D mesoscopic incompatibility in terms of the 
defect invariants for a countable set of dislocations. 

Theorem 2.4 (main result) For a countable set of parallel defect lines L and under 
assumptions and and definition incompatibility as defined by equations 

Sl.^) and fl.3\) is the vectorial first order distribution 

Vl = SkzV*z + SkK^K, (2.54) 
where its out- of -plane component is 

r,: = ("^'^''^lW + {Bf^ + ep^i&f - ^o/3)f^:''') (2.55) 

while its in-plane components are 

r?: = (2.56) 

and where xg £ ^ is a selected reference point. 

Proof. With use of theorems 12.21 and 12.31 the L^{Q) strain decomposes as: 

^mn ^mn ~^ ^mn^ (2-5T) 

where is compatible on Q, and where 

Emn ~ -E-mn (2.58) 

is smooth away from the defect set 'C = [J i'*'. 

iei 

Now, the local strain assumption of Part A is satisfied by each individual Emn ■ 
Then, the Frank and Burgers vectors of each isolated defect line L''-* are defined ac- 
cording to theorem II. II and definition ll.4l In a next step, the strain contributions Emn 
associated with these isolated defect lines are removed from the decomposition of the 
strain provided by (|2.8|l . (|2.42p and (|2.43p . as allowed by the absolute convergence 
1)2. 41[) of all their weight series. This operation defines a strain reminder whose defect 
lines are the accumulation lines (or the so-called derived set) of the initial defective set 
C, and the extraction procedure of isolated defect lines is repeated on this derived set, 
and then repeated again by transfinite recursion until the Frank and Burgers vectors 
are defined for each line of £ (and not only for the isolated lines). Finally, applying 
the main theorem of Part A to each L''\ summing the results on T C N and recalling 
that both and E'^n provide vanishing contributions to incompatibility, the proof 
is achieved. □ 

The remaining of this section consists in a reformulation of the above result in 
terms of the defect density tensors 0*^, A*j,, q*^, and which, in 2D, simplify in 
the 9J,A^,aJ and K*j tensors defined by (|1.10[) - (|1.16[) . However, the sums are now 
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performed on a countable ensemble £ of rectilinear dislocations, and and are 
Radon measures in view of the inequalities 

^|f7:«|<(X), ^||Bf'|l<~, (2.59) 
iei' iei' 

where I' refers to any bounded subset i £ X'} of Zo 

Theorem 2.5 For a countable set C of parallel defect lines and under assumptions 
\1.1\ and \l.S\ and defimtion \l.Sl the mesoscopic strain incompatibility writes as 

Vk = 3^kO*z + ea,Bdanli3, (2.60) 

or equivalently as 

Vk = S^kOl + ekaida^li. (2.61) 

Proof. Consider any straight dislocation L''^ £ C passing by f(^) = {x)^\zo). From 
theorem 12.41 the associated incompatibility writes as 

vf = + e^^ (^Bf> + ep^{xf - aro.s)^:*')) 3^5^,,,) 

+<5fe«ie««B:^'^ac<5i,„. (2.62) 

Taking into account (|l.lUp . p.ll|l . p.l5|l and p.l4|l for a single line, and the relation 

da {{Xf} - Xofl)Sj^(i)) = da (^{x'j^^ - ajQfl )(5^(i) ) = (£^'^ - X0|3)^a5J^^^) , 

it results that t)^'*' can be written in the formulations (|2.60p or H2.61[l . and the result 
follows after summation over X. □ 

3 Displacement and rotation fields in 2D incom- 
patible elasticity at mesoscopic scale 

3.1 Position of the problem 

The principal objective of the present work is to pave the way for a mathematically 
rigourous treatment of dislocation homogenization from meso- to macro-scale (Van 
Goethem and Dupret 2009b). To this end, this section elucidates the link between 
incompatibility (expressed as a combination of the concentrated defect densities) and 
the multivalued displacement and rotation fields. For the sake of completeness, most 
results are expressed using complete 3-D tensor components (with latin indices, cf. 
Part A) under the hypothesis of a 2D strain field (whose components do not depend 
on z). 

On the one hand, the mesoscopic fields O^. — Olf.,Al — A*j,,7;J = 77*^, = r/l^ and 
the contortion concentrated distributions on the defect lines which, as shown in 

Part A and SJ21 provide all the information on the dislocation and disclination densities 
and the strain incompatibility. 

On the other hand, by theorem 11.11 the rotation field is a multivalued function of 
index 1 obtained on flc by line integration of di"^u}l — diLul = ekpqdpSqi (cf. (|1.4[l ) 
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while the displacement field is a multivalued function of index-2 obtained on flc 
by recursive line integration of dj'''^dl'''^ul. = flj^'' {S^i +ujIi). The Frank tensor was 
introduced as a distributional field dicol defined on the entire fl, which coincides with 
dl'-'ujl on every defect-free region and to which additional distributional terms are 
added to let this tensor be related to the strain gradient by (|1.4|l . 

Direct analysis shows that displacement is not the most appropriate vector field 
to describe the dislocations since the Burgers field defined from p.5|) by line inte- 
gration on flc of i9;'^'&fc — dibX exhibits more useful properties. In particular, hi is a 
multivalued field of index 1 (compared to the less tractable index-2 multivaluedness 
of displacement) and the properties of the Frank tensor and vector and loI, and 
of the Burgers tensor and vector dib% and bj, show a clear analogy. 

Both the Burgers and the rotation field are defined by means of a Riemann foliation 
F (cf. Part A) in the sense that mappings of the following kind exist: 

fir ^ F "-H"' 

where F := {{x, #C) for every x G Qc and every curve C joining xo to x, with #C 
the equivalence class of all curves homotopic to C in flc} while V is the projection of 
F onto flc, in such a way that P(x-, #C) — x. Accordingly, the relationships between 
the multivalued fields cj^ and 6^ (defined on F together with the projection V), and 
the distributional fields diul. and dibl. (defined on Q) are very similar. 

Careful analysis however reveals an apparent contradiction between the expected 
meanings of the Frank and Burgers tensors and their mathematical properties. Theo- 
rem l2.4l first shows that in the absence of disclinations (Q,*J'''^ — 0) but in the presence of 
dislocations, the incompatibility rjl does not vanish (this resulting from non-vanishing 
coefficients multiplying the Dirac mass derivatives (i9a(5^(i) ) in (|2.55p and (|2.56p ). Since 
by (|1.8|) incompatibility is the curl of the Frank tensor, the latter is not curl-free and 
surprinsingly cannot be the distributional gradient of a single-valued rotation field oj^. 
in the absence of rotational defects (a situation where uj^. is expected to exist in the 
entire domain Q and not only on flc)- 

Secondly, the link between the Burgers tensor and vector has also to be clarified, 
but the situation is slightly difi'erent since no pure disclinations exist while there are 
pure dislocations. Indeed, it should be noticed that the Burgers vectors of the defect 
lines depend on the reference point xo in case their Frank tensor does not vanish. More 
precisely, when xq is changed to x'q, the following transformation rule applies on each 
defect line Z/*-'' in the general 3-D case: 

whereas the Frank vector flf, will remain invariant together with its scalar product 
with the Burgers vector. So, if in 2D the Frank vector of a given defect line L''-* does 
not vanish (f2z''' 7^ 0), an appropriate change of xq can always generate arbitrary 
values of the edge Burgers vector components Ba^^-' for this particular line. 

Suppose now xo can be selected in such a way that all the Burgers vectors -B^'''' 
vanish while the Frank vectors do not (fiz^*' ^ 0). Then by theorem 12.41 the incom- 
patibility r]l does not vanish, and hence by (|1.9p the Burgers tensor is not curl-free 
and cannot be the distributional gradient of a single-valued Burgers field fej in the 
absence of translational defects for this particular reference point xq. 
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3.2 Completed Frank and Burgers tensors 

In order to resolve the problem posed in i|3B.ll . the tensors djUJk and dibl are completed 
by appropriate concentrated effects within the defect lines, without however modifying 
their relationship with the multiple- valued Burgers and rotation fields on Qc- These 
tensors are called the completed Frank and Burgers tensors. 

Definition 3.1 

COMPLETED FRANK TENSOR SjUJ^ := djul-nlj, (3.1) 

COMPLETED BURGERS TENSOR Sjfefe ~ £kj + (^kpq{Xp - Xop)5j(^g ■ (3.2) 

The following theorems justify the introduction of the completed Frank and Burg- 
ers tensors. 

Theorem 3.1 In 2D, the 2nd-order-tensor distribution QjCdl verifies the relation: 

DISCLINATION DENSITY 6*fc = eujdiSjUjl. (3.3) 

Proof. This statement is a mere consequence of (|2.60|l and the relation H1.8|l expressing 
incompatibility as the curl of the Frank tensor. □ 

Theorem 3.2 In 2D, the 2nd-order-tensor distribution BjbX verifies the relation: 

DISLOCATION DENSITY A*fe = eujdiSjbl. (3.4) 

Proof. This statement directly follows from (fOll . (|3T|l and ([S^J- □ 

From the above analysis, it results that the curls of the completed Frank or Burgers 
tensors vanish in Q in the absence of rotational or translational line defects and that 
in these respective cases these tensors are equal to the gradients of existing single- 
valued rotation or Burgers vector fields. It is worth noting that the same concentrated 
contortion term is substracted from the Frank tensor in (|3.ip and (|3.2p in order 
to provide the completed Frank and Burgers tensors. 

3.3 Integral relations and Stokes' theorem 

Returning to 2D tensor notations, the following result restates the main theorem 12.41 
or 12.51 in terms of the completed Frank and Burgers tensors (|3.1fl and (|3.2|l and the 
associated multivalued rotation and Burgers vector fields. 

Theorem 3.3 Under assumvtions \l.l\ and \l.'A the mesoscopic strain incompatibility 
for a countable set C of rectilinear dislocations writes as 




with in Q,c , the index-1 multivalued fields loI^ and 6^ given by 




19 



Proof. This proposition directly results from the introduction of p. If) and (|3.2|) in 
the main theorem. □ 



Line integration of the completed Frank and Burgers tensors in flc therefore pro- 
vides the rotation and Burgers fields. When this integration is carried out on a loop 
enclosing a corresponding 2D area, the dislocation and disclination densities can them- 
selves be integrated on this area. 

Theorem 3.4 (Stokes' theorem for the completed defect tensors) Consider un- 
der assumptions \T71\ and \l,^ a countable set of dislocations and/or disclinations and 
a 2D open set S C Qzg perpendicular to the defect lines and bounded by the counter 
clockwise-oriented Jordan curve C C flc which encloses once the defect subset Lc '■= 
{L''^ L'*^ n S / 0, i el}. Then the following equalities hold: 

/ 9,3Wfcdx-/3 = / eapda^pujldS = / QldS = E f2**''(52fc, (3.5) 
Jc Js Js lW^Cc 

[ dpbldxfj = f e^pd^^pbldS= f AldS^ BfK (3.6) 

•'c Js Js LMeCc 

Proof. Since 1)2. 59p results from the assumptions, the dislocation and disclination 
densities are Radon measures on Q,~g and hence can be integrated on S. Then (|3.5[l 
and ((321) directly result from ((331), iflTOl) . ((1161) and ((LTT]) . □ 

Remark 3.1 The vector d/iLj* does not verify Stokes' theorem, neither in the classical 
sense, since eapdaSpi^t singular at x^^\ nor in a measure theoretical sense, since 
eapdaQpu}*^ is not a measure but a first-order distribution. The same remark can be 
made about the Burgers tensor. Nonetheless, as often observed in the literature, even 
in an inappropriate context the formal use of Stokes ' theorem here gives a correct final 
result. 



4 Concluding remarks 

This paper is part of a work devoted to the development of a mathematical theory to 
analyse dislocated single crystals at the meso-scale by combining distributions with 
multiple-valued kinematic fields. The distributions are concentrated along the defect 
lines which in turn form the branching lines of the multivalued fields. From this anal- 
ysis, a basic theorem relating the incompatibility tensor to the Burgers and Frank 
vectors of the dislocations and disclinations has been established in the case of count- 
ably many defect lines, under precise hypotheses on the distributional elastic strain 
gradient (via the Frank tensor). Quite surprisingly the sums of the norms of the 
Burgers and Frank vectors of the defect lines - which can be derived from the elastic 
strain - are required to be locally bounded to obtain the proof, thereby providing a 
fundamental defect norm for a further homogenization of the medium properties from 
meso- to macro-scale. This latter problem is addressed in Van Goethem & Dupret 
(2009b). 

Moreover, in addition to the elastic strain, two key objective internal fields (the 
completed Frank and Burgers tensors) have been identified to represent the medium 
defective state independently of the selection of the reference configuration. While the 
curls of these two first-order distributional tensors are precisely the disclination and 
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dislocation densities, their recursive line integration in the defect-free region provides 
the multiple valued rotation and displacement fields. 

After homogenization from meso- to macro-scale, no concentrated effects will re- 
main present anymore in the macroscopic model, which will consist of a set of evolution 
PDE's governing the tensorial defect densities in the framework of elasto- or visco- 
plasticity (cf e.g. Kratochvil & Dillon 1969). More precisely, the thermo-mechanical 
macroscopic model will govern the homogenized elastic strain and completed Frank 
and Burgers second-order tensors. Let us also mention that the non-vanishing meso- 
scopic elastic strain incompatibility will generate a macroscopic plastic strain which 
cannot be defined independently of the choice of the reference configuration. This 
property simply shows to be a reminiscence of the mesoscopic displacement and rota- 
tion multivaluedness. 
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